Forest conversion into pasture selects dung beetle traits at different biological scales depending on species pool composition

Abstract The conversion of forests into open areas has large effects on the diversity and structure of native communities. The intensity of these effects may vary between regions, depending on the existence of native species adapted to open habitats in the regional pool or the time since habitat change. We assess the differences in species richness and functional diversity of dung beetle communities (Coleoptera: Scarabaeinae) between native forests and novel pasturelands of the Atlantic Forest and the Cerrado, two biomes with contrasting histories of human occupation in Brazil. We conducted standardized surveys in seven forest fragments and adjacent pastures in each region and measured 14 traits in individuals collected in each type of habitat at each particular site. We calculated functional richness, functional evenness, functional divergence, and community‐weighted mean of traits for each area, and analyzed individual variation through nested variance decomposition and Trait Statistics. Communities were richer and more numerous at the Cerrado. We did not find any consistent relationship between functional diversity and forest conversion beyond the changes in species diversity. Although landscape changes were more recent at the Cerrado, the colonization of the new habitat by native species already adapted to open habitats lessens the functional loss in this biome. This indicates that habitat change's effects on trait diversity depend on the regional species pool rather than on time since land conversion. Forest conversion effects were primarily due to internal filtering. The effects of external filtering only appear at the intraspecific variance level, with contrasting differences between the Cerrado, where traits related to relocation behavior and size are selected, and the Atlantic Forest, where selection operates for traits related to relocation behavior and flight. These results evidence the importance of considering individual variance to address the responses of dung beetle communities to forest conversion.


| INTRODUC TI ON
Forest conversion is a major threat to biodiversity in tropical landscapes (Newbold et al., 2015). The conversion to open areas has large effects on native communities, through changes in habitat structure, the exclusion of native species, and the facilitation of invasions. Such replacement of native species by aliens may affect ecosystem functioning and decrease the effectiveness of the community in utilizing resources and resisting other disturbances (Harrison et al., 2014).
Trait-based functional diversity can provide means to assess these effects on the biodiversity-ecosystem functioning relationship (Flynn et al., 2011).
Functional diversity can be measured from the range of variation in the functional traits of the species present in the community, assuming that ecological functioning can be indirectly assessed through the diversity of traits with functional meaning (Dıáz & Cabido, 2001). Within this framework, a functional trait is any measurable characteristic of the individual (morphological, biochemical, phenological, physiological, and behavioral) that affects either its fitness, the fitness of other individuals of the same or different species, or other abiotic ecosystem processes (Violle et al., 2007).
Traits are used under the assumption that these characteristics provide information on the ability of individuals to perform particular functions and/or respond to biological interactions, thus providing a good proxy for ecological functionality. Therefore, by measuring different aspects of trait variation, different indices of functional diversity are thought to account for different aspects of functioning (Mason et al., 2005): Functional richness measures the functional (i.e., trait) space occupied by the species in the community; functional evenness does so for the regularity in the use of this space; and functional divergence accounts for how the differences in the distribution of the species in the trait space, which may contribute to better use of resources.
The use of traits in functional ecology is less common-though increasing-for animals than for plants and has been mostly focused on the study of assembly processes (Moretti et al., 2017). One common approach to conceptualize the assembly process is community filtering, where a series of filters determine which species are able to colonize the focal habitat fragment or locality and survive there. These filters are commonly divided into two categories: environmental (i.e., abiotic conditions) and biological (including competition, facilitation or density-dependent processes). Although this approach widely used, using it poses some challenges, as some biological filters may affect environmental filters, and vice versa. This led to the proposal of external and internal filtering processes (Violle et al., 2012). External filtering would select species from the pool on a scale larger than the community, through environmental or biological factors such as large-scale climatic gradients, or predator pressure along the landscape. Whereas internal filtering would encompass processes occurring locally within the studied community, like density-dependent processes or microclimatic heterogeneity.
These two types of filters are relative to the spatial scale of the community of interest, so this approach helps overcoming the complex interpretration of traditional filters, whose effects are often impossible to separate one from another (Violle et al., 2012). Nonetheless, this approach allows evaluating the scale at which the largest effects on the community are occurring through intraspecific trait variation.
Using only mean trait values per species does not allow assessing, for example, the effects on individuals who have trait values around the optimal mean, and can be benefited from density-dependent processes such as competition, increasing their fitness in the community. In this case, internal filters can increase variability by reducing the competitive pressure on these individuals, affecting the distribution of trait values around the mean optimal value selected by the external filter (Turcotte & Levine, 2016).
It is important to highlight the reduced number of experiments assessing trait functionality in animal functional ecology (see Noriega et al., 2018 for insects). Dung beetles (Coleoptera: Scarabaeinae) are to some extent an exception to this, being one of the few groups where several of these experiments had been carried out (deCastro-Arrazola et al., 2020;Emlen et al., 2005;Macagno et al., 2016;Nervo et al., 2014;Slade et al., 2007;but see deCastro-Arrazola et al., 2023).
Indeed, dung beetles can inform about the processes involved in the responses to forest conversion. They are a good model for these studies because they present rapid responses to ecological changes and are easy to collect (Gardner et al., 2008;Nichols et al., 2007).
Dung beetles are well-known for their feeding on mammal feces and the behavior of making and rolling dung balls shown by some of them (Halffter & Matthews, 1966). Their most iconic function is dung removal, but they also provide other functions such as parasite and fly control, soil bioturbation, contribute to diminishing CO 2 emission in pastures, incorporate NO 3 in the soil, and act as secondary dispersal of seeds and enhance plant growth (deCastro-Arrazola et al., 2023;Nichols et al., 2008;Slade et al., 2016). Their distinct feeding behaviors provide a classification in guilds that provides a rapid approach to their functional diversity (Doube, 1990; see also Pessôa et al., 2017). They can be classified as Rollers that make a dung ball and roll away; Tunnelers that burrow the dung; and Dwellers that K E Y W O R D S Atlantic Forest, Cerrado, community structure, functional diversity, habitat structure, land transformation, Scarabaeinae dung beetles, trait selection

T A X O N O M Y C L A S S I F I C A T I O N
Biodiversity ecology, Community ecology, Entomology, Functional ecology live directly in the dung (Bornemissza, 1969; see also Tonelli, 2021).
This knowledge of their natural history may help to better interpret the patterns observed in nature.
In the case of Neotropical dung beetles, the conversion of forest into pasture affects community structure by diminishing their richness and increasing the dominance of a few species (Nichols et al., 2007;Sánchez-de-Jesús et al., 2016). In functional terms, forest conversion affects dung beetle food relocation behavior, body size and daily activity (i.e., diurnal, nocturnal, or crepuscular) (Nichols et al., 2013), as well as their effects on ecosystem service provision (Noriega, March-Salas, et al., 2021). Although the effects of land-use change on the spatial and temporal dynamics of Neotropical dung beetle communities are relatively well-known (Dale et al., 1994;Gardner et al., 2008;Klein, 1989;Korasaki et al., 2013;Lopes et al., 2011;, there is a need for a better understanding of their responses considering their evolution in more forested areas, in comparison with Afrotropical and Palearctic regions. Functional diversity indices can inform about these responses, and dung beetle functional richness and divergence decrease as the impact of forest conversion increases (Barragán et al., 2011).
The Atlantic Forest and Cerrado biomes are both biodiversity hotspots, and their biotas are the result of distinct evolutionary histories and ecological processes that, arguably, have resulted in different regional pools of species. In general, dung beetle diversity in the tropics is greater in the forests than in open areas (Hanski & Cambefort, 1991; but see Silva et al., 2014). However, their ecological particularities create a conspicuous difference in the diversity of both biomes (Durães et al., 2005), since the natural landscapes of the Cerrado (aka. the Brazilian Savannah) host more natural open areas than those of the Atlantic Forest. Furthermore, the history of forest conversion in Brazil is spatially uneven (Leite et al., 2012). The Atlantic Forest was one of the first areas to be converted, mostly because it is situated in the coastal region, providing easy access to European settlers (Dean, 1997). Whereas the Cerrado was exploited more intensively in the expansion and internalization of the Brazilian population promoted by President Getulio Vargas in the 1950s (Oliveira & Marquis, 2002). Given the contrasting ecology and history of these two biomes, we expect that the differences in the functional adaptations evolved by dung beetles at each one of them would also affect their ability to colonize the novel open habitats.
We evaluate the effects of forest conversion into pasture on the functional diversity of dung beetle communities in the Atlantic Forest and the Cerrado. More specifically, we use data on community composition and trait measurements gathered from standardized surveys of forest fragments and pastures from seven landscapes within each biome, to evaluate whether forest conversion selects particular functional traits of dung beetles in each region through functional diversity indices and trait variations both between and within species. Therefore, we aim to answer the following questions: (1) How does forest conversion in the Atlantic Forest and the Cerrado affect the richness and functional structure of dung beetle communities? We expect that the effect in richness will be stronger than in the functional structure since dung beetles have a high functional redundancy; and that the time of conversion may decrease this effect (the Atlantic forest conversion event was about 100 years ago, whereas in the Cerrado, conversion happened roughly 40 years ago). (2) Is there a shift in the values of functional traits in the novel habitats created by the forest conversion? Here, we expect that traits related to dispersion or food reallocation will show larger values in pastures since the resource is more exposed, while in the forest traits related to maneuverability will be more important since the forest presents more obstacles during flight. (3) Which scale represents the variance found in the traits studied?
We expect that individual variations (i.e., intraspecific variance) and internal filters will be most important in the habitats with greater competition intensity.

| Study areas
This study was carried out in two different regions of Brazil: the Itajaí Valley (Santa Catarina) and the surroundings of Goiânia (Goiás) ( Figure 1). The Itajaí Valley is part of the Atlantic Forest biome, an evergreen tropical rainforest that has a constantly humid condition.
This biome comprises different vegetation types, such as seasonally semideciduous and deciduous forests, mixed ombrophilpus Araucaria forests, and ombrophilous dense forests (IBGE, 2012). In this study, we selected all fragments in ombrophilous dense forests, a formation characterized by large trees with dense crowns, which can reach 35 m in height, giving rise to a continuous canopy structure, and by a dense shrub understory, formed mainly by shrubs, herbs, and seedlings. There is also a wide variety of epiphytes, consisting of bromeliads, orchids, ferns, mosses, and lichens, resting on the branches of trees and shrubs (Vibrans & Sevegnani, 2013).
Goiânia region is located in the Brazilian Cerrado. This biome is subject to a regular and long drought period, which can last around 6 months. It is characterized by highly heterogeneous vegetation, composed of a continuum of areas of savanna, ranging from open grasslands, with no trees or shrubs ("campo limpo"), to forests (locally known as "cerradão") (IBGE, 2012). We selected all forest sites in fragments of "cerradão". This formation is characterized by the predominance of tree species, which can reach from 8 to 15 m or taller, giving rise to a continuous canopy. Many species are deciduous, so the crown cover can vary from 50% to 90% throughout the year. The understory is formed by small shrubs, herbs, and a few types of grass (Sano, 2008). Trees' crowns cast a considerable shadow, which makes the shrub and herbaceous layer smaller when compared to the other types of formations in the Cerrado.

| Dung beetle surveys
In each region, Cerrado and Atlantic Forest, we selected seven areas separated at least 1 km from each other. In each one of those areas, we conducted standardized surveys in two adjacent sites: one of forest and other of pasture. Dung beetle captures were made with baited pitfall traps consisting of 1-L pots with a solution of water, salt, and detergent. The baits were suspended above the trap with wire in a 50-mL plastic cup ( Figure A1). Three different types of baits were used: human feces, rotten liver, and cow dung. We placed three replicates of each type of bait, so in total, nine pitfall traps were placed in each sampling site, spaced 50 m apart along a transect, and separated at least 50 m from the edge of the fragment ( Figure A2).

The traps remained for 48 h in both habitats (forest and pasture).
We considered each pair of habitats as a sample unit. The surveys were conducted in the rainy seasons of 2016 and 2017. All collected beetles were identified by Fernando Z. Vaz-de-Mello (Universidade Federal de Mato Grosso) and deposited in the entomological collection of the Universidade Federal de Goiás. The dung beetle species pool in each region was obtained from the results of our surveys. Fragment size, shape, and conservation status may have some effects on dung beetle communities. We dealt with these undesired effects by surveying pasture and forest fragments in a pairwise manner. Also, we accounted for the particularities of each area by using it as a random factor in the models.

| Measuring dung beetle functional traits
We compiled information on a set of functional traits for each species and site based on measurements of the dung beetle individuals collected in the surveys. In total, we selected 15 traits ( Table 1) related to: dispersion (wing load, wing area/length ratio, and eye dorsal area; Byrne & Dacke, 2011;Dacke et al., 2013;Hongo, 2010); excavation (prosternum height, protibiae area, pronotum width, head length, and head width; Halffter & Matthews, 1966;Vilhelmsen et al., 2010); resource use (body size, measured as pronotum length + elytra length, and volume measured as length × pronotum width × prosternum height; Andresen, 2003;Emlen et al., 2005;Radtke & Williamson, 2005); food relocation (horizontal displacement and metatibia length; Halffter & Matthews, 1966); breeding behavior (nesting habit and nest shape-pear/ball; Halffter & Matthews, 1966), diel activity (Hernández, 2002); and specialization (i.e., food specificity; Falqueto et al., 2005). To obtain trait measurements, pictures of each individual were taken with a digital camera and using a stereoscope for smaller individuals, and the traits were measured in the software ImageJ (Rueden et al., 2017), using a graduated mm paper as a measure reference.
Food specificity was measured using Levin's index of niche breadth (Levins, 1968), based on the abundance of individuals of each species in traps with each type of bait of all traps placed in the same region, assuming that the wider the niche, the more generalist is the species.

| Functional diversity indices
Trait measurements were used to calculate three functional diversity indices: Functional richness (FRich), functional evenness (FEve), and functional divergence (FDiv) (Mason et al., 2005;Villéger et al., 2008). The indices were calculated for each habitat in each region, using the mean of the five individuals measured from that habitat. FRich measures the functional space of a community and is calculated by the convex hull volume of all the traits of the species present in the community. FEve represents the regularity of abundance of the species in the functional space and is measured by using a minimum spanning tree based on trait similarity between species or individuals. FDiv represents the degree to which the distribution of species in the functional space maximizes the divergence of traits in the community and is calculated by measuring the distance of the species to the centroid of the functional space. All functional indices were calculated using the "FD" function (Laliberté & Legendre, 2010).

| Data analysis
2.5.1 | Effects of forest conversion on species richness and functional diversity Taxonomic differences in community structure between habitats (forest and pasture) and between regions (Cerrado and Atlantic Forest) were characterized through principal coordinate analysis (PCoA) sample ordinations. We used generalized mixed models to assess the effects of habitat on species richness and all functional indices, accounting for area differences by including this factor as a random effect, with the following model: We also calculated the standardized effect size (SES) to remove the effect of species richness in the functional indices. For this, we first created null models by randomizing the community matrices with the function "randomizeMatrix" from the Picante package (Kembel et al., 2010). We set the null model for the "independent swap" algorithm proposed by Gotelli et al. (2011), which maintains species occurrence frequency and sample species richness. The  The CWM combines species trait data with abundance to assess the functional composition of assemblages, weighting the mean trait value of all species in an assemblage by their relative abundances.

| Taxonomic and spatial scales influencing dung-beetle communities
To To calculate the magnitude of the differences between the ob- served T-statistic values and those coming from a random assembly of individuals, we estimated the SES (Gotelli & McCabe, 2002) as follows: where I obs is the T-statistic observed value and I sim and sim are, respectively, the mean value and standard deviation of the null models (n = 999 randomizations). The null models were simulated with randomization procedures for each T-statistic. We used the "Cati" package (Taudiere & Violle, 2016) for calculating T-statistics, SES, and generating the null models.

| Effects of forest conversion on species richness and functional diversity
In total, 2681 individuals were captured in our surveys: 2143 in the Cerrado and 538 in the Atlantic Forest, pertaining to 63 species from 18 genera ( Table 2, Tables S1 and S2).
In the Atlantic Forest, the most common species were Canthon We found a regional effect and a marginal habitat effect on dung beetle richness ( species abundance data evidence the differentiation in the species pools of both regions (Figure 2). But also, that species composition differs clearly between forests and pastures in the Cerrado, while these habitat differences are smaller for the Atlantic Forest, since sites from both types of habitats largely overlap in these PcoA Axes ( Figure 2). None of the functional indices were related to forest conversion, even when removing the effects of richness and abundance by calculating the SES (Figure 3 and Table 3). Nonetheless, when traits are analyzed individually, we found regional effects in all traits except for the eye dorsal area and volume (Table 4). Note: S stands for species richness, FRich for functional richness, FEve for functional evenness, and FDiv for functional divergence. Significant and nearly significant models are highlighted in bold.

TA B L E 3
Results of the linear mixed models for the effects of habitat and region (and their interaction) on species richness and functional diversity indices, and their standardized effect sizes (SES).

F I G U R E 2
Pcoa Axis for the dung beetles surveyed in forest and pasture habitats of the Goiânia region and the Itajaí Valley (placed at the Cerrado and Atlantic Forest biomes, respectively). Circles represent Goiânia region, and triangles represent Itajaí Valley. In red Forest and in green Pasture.

| Forest conversion and shift of functional traits in the novel habitats
The analyses on individual traits obtained by CWM analysis showed that the values of most of them differ between regions, with the Atlantic Forest (Table S3) presenting larger values and greater variance in the community-weighted mean for all continuous traits (Table 4; Figure 4). In contrast, habitat type only showed nearly significant effects on wing load.

| Taxonomic and spatial scales influencing dung-beetle communities
Habitat contributed little to the nested variance of traits, and the differences between species were the principal factor that promoted variance in both regions ( Figure 5). Trait Statistics partly corrobo-

| DISCUSS ION
Our results show that the conversion of forest to pasture affected mainly species composition both in the Atlantic Forest and the Cerrado. Strikingly, although the overall functional structure of the communities was apparently not affected by habitat changes, the decomposition of these effects by traits points significant changes between regions and habitats. This is especially true for intraspecific variation, as the variance of traits between habitats and regions comes from species differences, fostered by internal (i.e., within species) filtering in almost all traits, with a small contribution of external filtering processes that affected only a few traits. rare and the pasture is currently used only by a handful of generalist native Atlantic Forest species and invaders. These differences in the effect of, arguably, the same type of habitat filtering are the consequence of the ecological and evolutionary differences throughout the historical formation of the pool of the two regional communities (sensu Ricklefs, 2015).
In the Argentinian forests, the regional context reflected in different degrees of impact of forest conversion on dung beetle communities, where humid forests presented a higher impact on dung beetle functional diversity than dry forests (Guerra Alonso et al., 2022). This contrasts with our results, where these regional differences did not reflect directly on the functional diversity of dung beetle communities. When the differences due to richness are removed (i.e., by using the Standardized Effect Size; Gotelli & McCabe, 2002), functional shifts between both types of habitats are relatively small, in apparent contrast with the large functional  Andean, Amazonian and Caribbean forests, Noriega, March-Salas, et al., 2021). This generalized loss of functionality has been attributed to the smaller diversity of resources (Lumaret et al., 1992) and/or changes in microclimatic conditions (Gómez et al., 2018).
The contrastingly smaller losses of trait variations found in our analysis compared with these studies may be due to our use of  et al., 2022). Besides that, the reduced loss of functional diversity in some biomes may be due to a high functional redundancy in their pool of dung beetle species, which allows maintaining the functionality in each type of habitat despite the regional changes in species composition. This difference between species is observed in the effects on trait community weighted means. All traits presented only a regional effect (evidencing the differences in the functional solutions present in the pool of each biome), while only wing load presented a marginal difference between pasture and forest. All traits presented a higher CWM in the Atlantic Forest, due to greater habitat heterogeneity presented in this habitat.
The lack of difference in the CWM of the traits contrasts Guerra Alonso et al. (2022) findings that traits related to size and food relocation presented clear differences between forest and pasture, mainly due to forest energy restrictions to telecoprids and the consequence dominance of this habitat by paracoprids.
The partition of trait variance shows that interspecific variation had a greater contribution for the total variance of traits in both regions, corroborating the results found both in dung beetles (Griffiths et al., 2016) and other groups (De Bello et al., 2009;Messier et al., 2010). However, even though traits vary more between than within species, our results indicate significant intraspecific variation between habitats. In fact, our results show no effects of filters when ignoring intraspecific variation (the metric T PC∕PR ), while showing signal of external filtering for some traits when considering intraspecific trait variation (the metric T IC∕IR ). This result emphasizes the importance of considering intraspecific variation in community studies (MacArthur & Levins, 1967;Violle et al., 2012).
Eye Dorsal Area in both regions, and Metatibia and Head Length in the Itajaí Valley, presented a slightly higher intraspecific contribution than other traits, although still much smaller than interspecific variation. Eye dorsal area can be related to both flight ability, a period of daily activity, and the adaptation to different light conditions (Byrne & Dacke, 2011 the forest probably because of the highest number of rollers and bigger species than in the pasture. In this biome, we also found intraspecific variance in the prosternum height and protibia area, two related to excavation (deCastro-Arrazola et al., 2023;Halffter & Matthews, 1966). Indeed, soil texture and compaction affect the assembly of dung beetle communities (Davis, 1996). Therefore, the uneven compaction of the soil in the pasture may be selecting a larger interspecific variance in these traits, through the selection of individuals adapted to exploit soils both well-developed soils and those that have been compacted by cattle.
The greater promoter of individual variance in our data is internal filtering, which is consistent with other studies that use Tstatistics (Luo et al., 2016;Mungee & Athreya, 2021;Xavier Jordani et al., 2019;Zorger et al., 2019). The metric T IP∕IC was lower than expected by chance for almost all traits, suggesting little niche overlap, which can be promoted by local processes such as competition.
Indeed, those strong internal filtering effects were expected, since dung beetles present highly competitive communities due to the use of an ephemeral resource (sensu Atkinson & Shorrocks, 1981;Elton, 1966), which may be even stronger in the pasture considering that microclimatic conditions of the dung pat diminish the opportunity window of resource availability.
Several traits show a signal of external filtering, presenting opposing patterns in the two habitats. While the external filtering processes of the forest promoted overdispersion, increase in niche overlap, in the pasture they promoted clustering, niche packeting.
In the pasture, the fluctuation of heat and humidity may impose an important filter for selecting species and individuals with particular trait values. While in the forest the greater environmental stability promotes heterogeneity in the traits and the persistence of more strategies for resource utilization. At the Cerrado, forest habitats increase the individual variation of Body Length and Volume, and metatibia length, while in the pasture, the individual variation in those traits and pronotum width decreases. In the case of metatibia length, a trait related to the ability to roll dung balls (Halffter & Matthews, 1966;Hanski & Cambefort, 1991), this effect may be due to the lower presence of rollers in the forest (Krell et al., 2003). The dominance of tunnelers and dwellers in the forest may increase the individual variation in this trait, in contrast to the dominance of rollers in the pasture. Length, Volume, and Pronotum width represent different aspects of body size. Finally, individual variation in size may determine the amount of resources utilized for development (Emlen et al., 2005). Therefore, the greater variation in the forest may be a reflection of the uneven availability of resources in contrast with the  We thank Fernando Z. Vaz de Mello for dung beetle species identification.

CO N FLI C T O F I NTE R E S T S TATE M E NT
The authors declare no conflict of interest.

DATA AVA I L A B I L I T Y S TAT E M E N T
The data of the present work can be found as supporting information in the online version of this paper.

F I G U R E A 2
Design of the surveys of the dung beetle communities in two regions (placed at the Atlantic Forest and Cerrado biomes). In each region, seven areas were selected, and in each area, two habitats were sampled: one forest patch and one pasture adjacent to the forest. In each habitat (forest and pasture), nine pitfall traps with three different types of baits (human feces, rotten liver, and cow dung) were placed along a linear transect at 50-m intervals. In both transects, traps were placed at least 50 m from the habitat edge, to account for edge effects.

F I G U R E A 3
Dung beetle functional traits measured in five individuals per habitat per area. The traits were measured by ImageJ. 1. Dorsal eye area, 2. Head length, 3. Head width, 4. Pronotum length, 5. Pronotum width, 6. Elytra length, 7. Protibia area, 8. Metatibia length, 9. Prosternum height, and 10. Wing area. Body length was calculated summing pronotum length and elytra length. Wing load was calculated by the ratio of wing area by body size. And volume was calculated by multiplying body size, pronotum width, and prosternum height.